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ABSTRACT 


Empirically derived expressions describing the temperature-Structure 


A Te ; 
parameter, C_ , in terms of the stability parameter zZ/L or the Richard- 


it 
son number, Ri, over land were used to analyze data of mean humidity, 
temperature and wind speed and temperature fluctuations obtained by ship- 
board measurements over open ocean conditions. 

In general, there was little correlation between the spectral 
analyzed ene values and Z/L or Ri. It is recommended that a different 
approach be used to obtain ong estimates and that z/L or Ri be modified 


to account for the anomalous temperature behavior found over a marine 


environment. 
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im INTRODUCTION 


The advent of more complex and complete numerical models for analyz- 
ing and predicting large scale atmospheric motion has also led to better 
specifications of the boundary layer. The boundary layer is that region 
in the first kilometer over the sea which is defined principally by turbu- 
lent transfer of momentum and heat. Progress has been made in specifica- 
tions of mean surface wind fields and hydrostatic stability conditions. 
These are important parameters for predicting changes in the upper part 
of the ocean or for estimating boundary fluxes in medium range atmospheric 
prediction models. 

The Navy, however, also requires surface layer specifications with 
regard to their effect on the propagation of optical waves, and the 
Naval Postgraduate School has been tasked with providing the basis for 
prediction of the optical properties of the atmosphere for relevant atmos- 
pheric window wavelengths. The objective of such a program would be to 
relate the optical properties and the optically relevant meteorological 
parameters of the marine boundary layer to the observable bulk atmospheric 
parameters. 

The micrometeorological quantities which define the bulk meteorologi- 
cal quantities are time-averaged values of wind, temperature and relative 
humidity measured at multiple levels within the constant flux layer. 
Optically relevant parameters measured are turbulence induced temperature 
and humidity fluctuations, their frequency spectra, and the aerosol and 
fog concentrations. The optical properties measured are the modulation 


transfer function of the atmosphere, the image wander, its frequency 
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spectrum, the extinction due to absorption and aerosol/fog scattering, 
and the optical turbulence-structure function, ea related to scintilla- 
tion and defined in more detail in section II. It becomes of prime im- 
portance in the investigation of optical properties in a marine 
environment to compute easily oe Fortunately, a relationship exists 
between ei and the temperature-structure function, Chae: which is derived 
from micrometeorological considerations. 

A serious handicap in describing the marine environment in the past 
has been the lack of data and the fact that available expressions, which 
are empirical in nature, have only been validated over land. 

This study is an examination of data collected over a period of two 
years under open ocean conditions from the R/V ACANIA. The data are com- 
puted with expressions derived by Wyngaard et al (1971) and Friehe (1976). 
These expressions relate og to a stability parameter, either the 
Richardson number or z/L. Both the temperature-structure function and 


the stability parameters are defined later in this manuscript. 
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Ii, THEORETICAL BACKGROUND 


A. GENERAL 

Turbulence properties of interest in optical propagation are those 
which could be used to describe the intensity of fluctuations in the 
refractive index and that which could be used to describe the size or 
scale of the refractive index inhomogeneities. The former has bcen re- 
lated empirically to scintillation, beam spread and beam wander; and 
the latter along with the former to the image resolution which is im- 
portant with respect to seeing conditions. 

On the basis of the isotropic nature of small scale fluctuation, 
Only one parameter is necessary to describe the intensity of the refrac- 
tive index fluctuation over many scales. It is the refractive index- 


structure function parameter, previously mentioned, 


Ze 


S = [n(x) - SUERTE aaa, 


a) 


where n(x) and n(x+r) are refractive indices at two points ona line 
Oriented normal to the mean wind direction and separated by the distance 
r which is less than the outer scale, Loe (lower end of the inertial 
subrange) and greater than the inner scale, Xe (associated with the 
upper limit of the inertial subrange). The inertial subrange in this 
spectra is the wave number region where energy is transferred without 
dissipation from lower to higher wave numbers. See Figure l. 

The refractive index is determined primarily by density fluctuations 


’ oo 
and can, therefore, be related to temperature fluctuations. C. 1S 


1 
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Figure 1. A schematic drawing showing spectral transfer leading to 
existence of the inertial subrange (Lumley and Panofsky, 
1964). 


14 





2 
related to the temperature-structure function parameter, Cr , by the 


relation: 


jie (2) 


where os can be defined by the temperature variance spectrum, below. 
This turbulence parameter is the most important variable in determining 
the optical propagation characteristics of the atmosphere and will be 
discussed in the next section. 

The following expression for the one-dimensional spectral density 
representation of temperature fluctuation is derived from similarity 


theory and dimensional analysis: 


25/2 
1 — 
S, (x) Oa Cor is 


(3) 

As described earlier, turbulence is nearly synonymous with tempera- 
ture fluctuations; hence it is necessary to describe the mean thermal 
stratification. Parameters for this are the Richardson number defined 
by profile measurements or (z/L) defined by flux measurements. These 


concepts will also be expanded ina later section. 


Pee uCULATION OF THE TEMPERATURE-STRUCTURE PARAMETER 

There are several methods available to estimate a directly sand an 
directly. These have been partially introduced in previous discussions. 
The first approach is used to determine ae directly by the 2/3-law 
expression for the structure function. In a locally isotropic field, 


the temperature-structure function has the form 


C = [T(x) - nee 


ie (4) 


ES 





at separation r of the order of inertial subrange scales. 
A second approach is based on a one-dimensional temperature spec-~ 
Lreum, Sn (Kk), which by definition is the Fourier Transform of the corre- 


lation function with Separation r in the streamwise direction: 


a 5/3 ’ 
Ch = 4 S,,(k) k (5) 


where k is the streamwise component of wavenumber. Since velocity and 
temperature fluctuations are measured at a fixed point in the flow, the 
resultant spectra are realized at a temporal frequency, £. To obtain 

2 


Cr , the temporal (f) and space (k) scales are assumed to be related by 


Taylor's frozen-field hypothesis, i.e., 
k = 2nf/U (6) 


where U is the mean wind Speed at the measurement level. "Frozen" 
implies that the turbulence pattern remains unchanged as it Sweeps past 
the probe. 

The third approach involves measuring separately the rate of dissipa~ 
tion of turbulent kinetic energy (€) and the rate of dissipation of 
temperature variance (N) and applying the results to an expression 
Suggested by Corrsin (1951), as another form of the temperature spectrum, 


Viaiszey, 


Vee - 5/3 


a) BNE k (Ja) 


which yields the following expression for oe based on its definition in 
Eguation (5): 


2 


co = Bue.’ - 


(7b) 
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(8 is an empirical constant evaluated to be 0.25). From this last 
— 2 . . 
form, indirect estimates of on can be obtained based on mean condi-~ 
tions. This is because € and N are easily related to boundary fluxes 
and profiles, if steady and horizontally homogeneous conditions exist. 
zs ; 
Expressions which relate Cw to mean properties of the boundary laver, 
such as Ri or z/L, are desirable because the small scale measurements 
are impractical to obtain in most operational or tactical situations. 
2 or : 
Expressions relating Ch to € and N, Equation (7a), and similarity 
theory predictions for the dependence of € and N on momentum and heat 


fluxes were obtained from examinations by Wyngaard et al (1971): 





3 = FE (2/) (8) 
OM Z, 
N 
ee 2 (2/15) (39 
U4 2 
t | Ve 1 t = 3 1 t 
where = = (-u'w') ae wT pus L = “TU, /kg w'T' and £, (2/L) 


and £, (2/0) are empirically determined functions. 
' Direct substitution of Equations (8) and (9) into Equation (7a) 


yields an expression of the form 


gee 2 2-2/2 
eo t, 2 £ , (2/L) (10) 


where £(2/L) results from the combination of £, (2/1) and £,(z/L). 
Furthermore, Since zZ/L and Ri can be functionally related, a parallel 


dependence on Ri can be obtained, which is (Wyngaard et al) 


C 2 7A/ 


3°.- 2: 
a (dT/dZ) £, (Ri) (ii) 


i7 





This final expression provides the desired dependence of a on 


more readily measured parameters (z, dT/8Z, and Ri). 


Cy CONSIDERATION OF THE RICHARDSON NUMBER 
1. Derivation of the Richardson Number 

McIntosh and Thom (1972) provide a clear and concise explana- 
tion of the Richardson number which is the ratio of the buoyant kinetic 
energy (KE), to the natural kinetic energy of turbulence (K.E.) of 
eee particle in a turbulent boundary layer. 

Two reference levels (2) and Z5) are chosen within a turbulent 
boundary layer (shown in typical profiles of velocity, temperature and 
Gemeity 1n Figure 2), such that, Zo = 2) + 2, where & is the local mix- 
ing length. In other words, turbulent exchange of each property occurs 
through the same characteristic distance, cr mixing length, &. 

The total kinetic energy of a fluid particle of volume V within 
this boundary layer is made up of two independent components: one asso- 
Ciated with the horizontal progress of the particle as part of the mean 
flow, and the other associated with its own random movements as part of 
the forced turbulence present in the boundary layer. (The latter is 
related to the work required to raise the particle from Z4 ice) Z5 by 
buoyant action and to the kinetic energy associated with the purely 
vertical movements of the particle.) 

The instantaneous velocity of the particle exceeds the mean value 
Pyeeeaneamount U' given by U(z.) - U(2Z,)i 1,6. , CO aufirst aboroximacionm, 


Meee (00/02). McIntosh and Thom also assume that it is justifiable 


to equate U' and w', so that the vertical component of the 'natural' 
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kinetic energy of turbulence, (KE.) can be written 


jes ee... oe ns Z Z 
(K.E.) | = 5 OVw = 5 OVU ai PVR (dU/0Z) (12) 


Due to turbulent exchange, the fluid particle initially at Z4 
and with density 4 appropriate to that level, will find itself at 
level Z. and in surroundings of density Po. At level Zoe therefore, 
the particle is subject to a positive buoyancy force of magnitude 
g(p, - p,)V =F. Assuming that a linear temperature profile exists 
locally between z, and z, and that the buoyancy force acting on the 


M Z 


particle increases uniformly from zero at Zy OE a: Zor the work done 
(W) by buoyant action is equal to the product of the force =F and the 


distance (Z, - Z1): 


Ik 
Say g(P., = p,? (Z,, = Z,)V (1:35 


This amount of work goes directly to increase the upward 
kinetic energy of the fluid particle, equal to (KeE.) by an amount 


(K.E.)), to which W may be equated: 
(Keeeee= = a(p. - Pp.) (z,- z.)¥ (14) 
meee oo An A] 


30 
minally, on Substitution of & for (z., - Zi)s R ‘yz for (0, p,) 


and -£2 for <2, assuming an incompressible fluid, the equation becomes 


é (15) 


Haha 
S18 


i 
(K.Eo)) on aad >, OV 


The magnitude of the ratio (K.E.) /(K-E.) | indicates the rela- 


tive importance of free and forced convection in determining the 
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structure of a fluid boundary layer (shape of velocity profile, nature 
OL turbulence, Cen tne Negative Of this ratio, or the Richardson 
number Ri, iS a measure of the stability or instability in a region of 
boundary-layer flow. In lapse conditions Ri is negative, and in inver- 
Sion COnditions, it iS positive. After cancelling appropriate terms 


_g ol 7oZ 


(16) 
ee (9U/3z) - 


which is one of the forms used in this study. 
PeeeOtner Relationships for Richardson Number 

The Richardson number Ri can also be expressed in terms of 
(276) Or as a bulk aerodynamic formula. 

The stability of the surface layer regime determines flux in- 
tensities. Monin-Obukhov theory determines the scaling level, L, which 
1S proportional to the level where mechanical and thermal production of 
turbulent kinetic energy are equal. The ratio of the height of the 
measurement (z) to the Monin-Obukhov length (L) serves as a stability 


index (zZ/L). The Monin-Obukhov length is defined as 


TU," 
ys 
a 17 
qk tT, a) 
V 
aoe | hh T,)!hU+ OL.61 gg, T, Ge 7 O70, and ih = 0335 (Von Karman 


Vv 
constant). 


In unstable conditions the stability parameter z2/L is approxi- 
mately equal to the Richardson number Ri. The following relationships 


between the stability parameter and Ri have been proposed by Dyer and 


Zi 





Hacks (1970) for an unstable case and Webb (1970) for the stable case: 


Ri 
Z2/L = 7 Gens (stable) (18) 


Z/L = Ri (unstable) (19) 


Here @ iS an empirically derived constent equal to 0.5. 

A bulk aerodynamic formula is one which relates fluxes across 
the air/sea boundary to the wind speed at a level and the temperature 
difference between that level and the surface. The derivation of such 
a formula involves several assumptions regarding the stability condi- 
tions of the boundary layer and the turbulent processes within it. If 
valid, however, this type of formula is very useful for most practical 
needs. 

Friehe (1976) selected bulk aerodynamic formulae for empirical 
parameterizations of the surface fluxes of momentum, sensible heat and 
moisture in terms of average surface layer variables. He used the 
selected bulk aerodynamic formulae to derive an expression for the 
estimation of ei in terms of the relatively easily-measured quantities 
mean wind, U, temperature, T, and water vapor density, Q. This formu- 


lation by Friehe is presented in the following section. 


Deeeoasbol ING EXPRESSIONS AND RESULTS 
Wyngaard et al (1971) examined forms of nondimensional expressions 
: : ae Zz 
relating atmospheric stability to the structure parameter, Cn Damen] 


Peedateanalysis Of the height (z), the friction velocity (U,), the 


= 2 
mean potential temperature gradient (96/92) 7, €, N and Co produce 


1 : 
In this study virtual potential temperature gradients were used 
to account for humidity gradients over open ocean conditions. 


a2 





functional expressions dependent only on the Richardson number. Simi- 
larity Mrtined pals are elucidated by using the following Richardson 
number definition to indirectly obtain a comparable temperature-struc- 
ture parameter, 1.e€., 


g/T (30/92) 


Ri = = = 
(dU/dZ) 


(20) 
This relationship differs from Equation (16) in that the temperature 
parameter here iS a mean potential temperature. When expressed in 
the finite difference form, a simple approximation to dU/oZ, for 


example, is 


Fe ue ao) 


oz IZ i2Z Z. in(z,/Z,) 


where the computed number applies to a geometric mean height defined 


= WL: 
by Z. = (212.5) Ne 


Equations including the Richardson number functions (£1. f., and 
£,) are equations (8), (9), and (11). Equation (11) relates the tem- 
perature-structure parameter and stability parameter and allows for an 


empirical determination of the £, function at varying conditions of 


} 


Peemitty. Lf the function fy is defined, the temperature-structure 
parameter can be readily calculated from more easily measured para- 
meters, 00/8Z and Ri. In Figure 3 the predicted curve and the data 
obtained from the extensive 1968 AFCRL study of turbulence structure 
over a flat, unobstructed Kansas plain are presented. As will be shown 


in the section describing results, available marine data does not appear 


to agree with this predicted curve. 


as 
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The dimensionless temperature-structure parameter vs. the 
Richardson number (Wyngaard et al, 1971). 
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Friehe (1976) combined the semi-empirical eguation of Wyngaard et 
al with the bulk aerodynamic formulae to provide alternate equations 
for the estimation of ee in the atmospheric boundary layer above the 
sea. These equations relate en to measurable quantities of average 
wind speed and air/sea temperature and moisture differences. 

The expression derived by Wyngaard et al and considered by Friehe 


is Eguation (10) with £. defined as 


8 
=2/3 
f. aero (il = 7 (2/1) } / FOr 4S 0 
f. =eaeo(i + 2.75(2/1L) ) EOL 2/1, 7 0 
and z = height, m 
T,. = characteristic temperature scale, K 


The bulk aerodynamic formulae considered by Friehe were 


te =¢. U- (22) 

wo = Craeeeeets) = C. UAT (235) 

wa = eo U(Q. oe) = Cc. uAO (24) 
where U = average wind speed at height h, m/s 

Ts = average potential air temperature at height h, K 

2. = average water vapor density at height h, ee 

a. = average sea surface temperature 

2. = average water vapor density adjacent to the sea surface, 


usually obtained by assuming the air to be saturated 
1.e., at dew point et g/m> 


i = drag coefficient, from experiment 
C. = sensible heat flux coefficient, from experiment 
= moisture flux coefficient, from experiment 
h = reference height, usually 10m above the average sea 
surface 
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w = vertical velocity fluctuation, m/s 
, S 
Gq = moisture fluctuation, g/m 
6 = potential temperature fluctuation, K 


The resulting expressions derived by Friehe and adapted for this 


study were the following: 


For -20 < UAT < 25 (mK/s): 


ee - 
Ske x LO Cite eyo) We unstable 
2 2 22/3 
| 
2 -3 
car) Bele xO (in=“O- 635%) stable 
where X = ae tp 0 2a R, AO/7 AT) Reem 
U 
For UAT > 25 (mK/s): 
pa 2/3 
&: Ze 
E Be eoeeion «(1 + 2.60x1) > 
2 
(AT) 
Za 
where X' = Sean! +eOets2 Ry AO7AT) 
U 


(25) 


(26) 


R., and R, are the ratios of the coefficients used in making adjust- 


3 4 


ments to the levels used for this presentation. Results obtained by 


Friehe (1976) appear in Figure 4. 
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Pi eee ear eERIMEN 


A. THE PLATFORM AND LOCATION 

Observations were obtained aboard the R/V ACANIA anchored off 
Monterey, California in Monterey Bay. Measurements were made at four 
levels on two masts spatially separated on the forward deck of the 
ship. An illustration of the sensor locations of these masts is shown 
in Figure 5. The vane and probe arrangement appears in Figure 6. The 
vane maintained a one-dimensional profile of the wind during measure- 
ments. 

The location of Monterey Bay provided an ideal site for the experi- 
mental program as shown in Figure 7. Open ocean differs from land in 
Mmemwerrects Of wave action on turbulence, in the nature of the aerosols 
and fog, and to a lesser degree but perhaps significantly on the humid- 
ity fluctuations. These conditions can, of course, be best obtained 
far at sea. The cost of such activities makes it desirable to work 
near land. Pt. Pinos and Monterey Bay provide very nearly the ideal 
eieeuvation. 

Pt. Pinos projects northwest from the mainland toward the prevailing 
northwest wind. Even under storm wind conditions, the wind comes from 
the southwest, still bringing sea air toward shore. The ocean depth sur- 
rounding Pt. Pinos is shallow enough that the ACANIA can anchor at any 
range up to about 15 km to the north and northeast. 

The effects of land in producing acceleration of the airflow are 
very small in the Pt. Pinos-Monterey environment. The shoreline south- 


east from Pt. Pinos is parallel to the prevailing wind. It rises 
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Figure 5. Mounting arrangements aboard the ACANIA. 
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Figure 7. Experimental site in Monterey Bay. 
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slowly along the water's edge. The eastern edge of Monterey Bay con- 
sists of stabilized sand dunes. The sand beach shore is the mouth of 
the Salinas Valley, about 20 miles wide at the shore line. This valley 
extends approximately 100 miles to the southeast. In daylight hours, 
it adds a convective airflow to the southeast. The conventional day- 
night sea breeze reversal often roughly cancels the prevailing airflow 
pemeame tO produce a calm period rather than a reversal, so that a 


relatively small amount of land air is transported out over the water. 


B. INSTRUMENTATION 
1. Mean Measurement Systems 

The mean wind measurements were made with a Thornthwaite Asso- 
ciates cup anemometer wind profile register system, model number 104. 
In operation the shaft of a three cup anemometer unit serves as the 
shutter between a light source and a photocell for each Gevolution. 
The cups are plastic cones reinforced with aluminum frames. They are 
attached to the rotating shaft by stainless steel tubing spaced at 
120 degree intervals about the shaft as shown in Figure 8. The three 
cup assembly sets along with the other sensors were positioned at roux 
levels on the bow with electrical leads to the after deck house labora- 
tory. The sets have the characteristics of low starting speeds with 
a small amount of internal friction which aids in checking inertial 
overshoot. 

Temperature sensitive quartz crystal probes, (Hewlett Packard 
model HP-2850) were used to measure mean temperatures at each level. 
RF signals from the crystal probes and from a reference oscillator were 


masa in the HP-2801A readout unit to produce a beat frequency whose 
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Figure 8. C. W. Thornthwaite Anemometer Cups. 
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Signature can be analyzed to within 0.001 degrees centigrade per hertz. 
Each sensor simultaneously received pre-experiment calibration against 

a platinum resistance wire thermometer in a temperature controlled cir- 
culating water bath over the expected temperature range. The accuracy 
in achieving a 0.005 degree centrigrade correction factor was a constant 
for each probe. A 3.7 meter flexible coaxial cable is permanently 
attached to the sensor head and the mast mounted probes are Helsed 19) 
an aspirated shelter as depicted in Figure 9. Temperature values were 
automatically recorded on a printer tape. 

The Hygrodynamics Digital I hygrometer indicator and the Dunmore- 
lithium chloride sensor were used to measure the relative humidity. The 
Dunmore-type lithium chloride sensor is also pictured in Figure 9. The 
equipment operates on the basis of resistance change within an electro- 
lytic solution generating a reference voltage variance which is propor- 
tional to the relative humidity change. Automatic temperature compensation 
in the instruments meets the following specifications for relative humidity: 

mae Leldative humidityspelow 90% 
4 4% relative humidity above 90% 
Calibration of the sensors was accomplished by comparative methods using 
Saturated salt solutions in an enclosed chamber. Relative humidity was 
observed and recorded similar to that of temperature as a printed output. 
2. Temperature Fluctuation Sensor 

The temperature fluctuations were measured using a bridge 
developed by personnel at GTE Sylvania, the GTE Sylvania Model 140. 

The system was slightly modified for use in this study. 

This system is a small temperature difference sensing device 


originally designed to be carried aloft by a balloon and equipped to 
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process and transmit, by pulse-rate modulation on a radio-frequency 
carrier, this information to a ground-based receiving station for 
demodulation and recording. 

The baseband portion of this system is basically a balanced 
Wheatstone bridge excited by a 3 kHz signal with a synchronous detector 
on the output. Segments of a very small diameter platinum wire serve 
as temperature sensors in opposite arms of the bridge. In the single 
wire mode, one wire is replaced by fixed resistance. The resistance- 
temperature coefficients result in an output from the bridge which is 
proportional to the temperature difference between the two probes or 
a temperature change of a single wire in the single wire mode. 

The sensor wire is 0.5 cm long and 2.5 x 10° m in diameter. 
This extremely small mass allows a response to temperature variations 
of up to 1 kHz while electronic amplification allows temperature differ- 


ences as small as 0.004 degrees C to be observed. 


C. ANALYSIS PROCEDURES 

This particular study was divided into two phases: 1) individual 
mean wind U, mean temperature T, and mean humidity gq values were plotted 
Guelogarithmic scale. Best fit lines were drawn to the points from 
mmenmedata were picked off and applied to the expressions developed by 
Wyngaard et al; and 2) values of U, T, and gq from individual levels were 
applied to the bulk aerodynamic expressions examined by Friehe. 

Migactuating temperature T’ and profiles of U, T and q were studied 
for eight time periods. Four levels of data were desired, but signals 
were not always available for all levels due to various component 
failures. Table I summarizes the periods examined and what data were 


available for each period. 
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SUMMARY OF DATA PERIODS 


NUMBER 
Osa 


LEVELS 


NUMBER 
OF 
PROFILES 


(T,U,q) 


12 


24 
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24 


2g 





Mean wind measurements over ten minute intervals were recorded in 
a log book. The recordings were the number of cup revolutions for the 
interval as recorded by electronic counters. Temperature readings were 
printed on a paper tape from a monitoring digital readout at approxi- 
mately one second intervals. Relative humidity measurements were re- 
corded the same as the temperatures. All information was later coded 
and punched on computer cards for further analysis. 

The data were then screened for gross errors or inconsistencies 
due to instrumentation malfunctions. The criterion at this point for 
retaining or discarding profiles or data periods depended on whether 
en values could be computed from temperature spectra for the time 
periods involved. 

The computer analysis of both phases consisted of programming the 
governing equations uSing either actual data or data picked off the 
profiles. The NPS Fortran program PLOTP was used to plot the predicted 
curves for phases one and two as well as the results from the data. 

fee Profile Editing 

The data from wi; T, and gq measurements were plotted on 4-cycle 
semi-logarithmic paper. Since U, T, and gq are parameters which vary 
logarithmically with height, a best fit straight line was drawn to the 
data points. In general, the procedures were subjective. In many in- 
stances several different slopes could be obtained from just one graph, 
so a useful criterion was not to give a Single data point too much 
influence in determining the straight line. Consequently, the line 
drawn represented a most probable position between data points as 


illustrated in Figure 10. 
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Difficulty arose whenever anomalous points appeared in a graph 
as in Figure ll where the temperature at the third level was obviously 
inconsistent with the second and fourth levels. In a case such as this, 
the point was merely ignored. 

The time period of 21 November 1974 (1817-2130) was particularly 
@ifficult to analyze due to the apparent discrepancy of the temperature 
and humidity sensors at the fourth level from the readings of the other 
levels. This discrepancy is evident in Figure 12. These sensors were 
found in later tests to be out of calibration, and the decision was 
made to draw the straight line to fit the three lowest levels. 

In some, but not in all, cases the sea surface temperature was 
helpful in determining the slope. However, if the sea surface tempera- 
ture was significantly different from that expected on the basis of the 
temperatures in an atmospheric profile, it could not be used since the 
profile temperatures differed in most cases by only hundredths of a 
degree. 

2. Analysis Procedures to Obtain oe 

Variance spectra were interpreted to estimate het The spectra 
were obtained using an analog spectral analyzer. Procedures for convert- 
ing temperature spectral values, obtained with the analog spectral 
analyzer, to engineering units and then for obtaining the turbulence 
parameter eae from the temperature spectral estimates are described 
Pees Section. 

a. Spectral Plot Scaling 

A necessary procedure was to scale the spectral plots in 
order to relate RMS input voltages to power spectral densities; (vari- 


ance per unit frequency). 
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Figure 12. Profiles of 21 Nov 74 (1817). Solid black line represents 
the profile used. 
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To obtain power spectral density levels, corresponding to 
RMS voltage inputs, calibrated scale charts had to be constructed. The 
charts were constructed as follows. Amplitude scaling was accomplished 
by using an externally generated signal whose RMS value was determined 
by a RMS meter. Using a O-dB (.lv) input gain and a spectral gain of 
O-dB (x 1), a signal with amplitude equal to 0.1 volts RMS was produced 
on a screen display as a spike near the selected frequency. This spike 
was then plotted with a xX-Y¥Y plotter. Successively, the input gain was 
stepped down to attenuate the amplitude of the input Signal by 10 dB 
increments, and plots of the height of each resulting spectral spike or 
amplitude was then added to the xX-Y plot, at different frequencies. 
These plots represented a graduation of RMS input from a minimum of 
0.001 volts RMS to a maximum of 0.1 volts RMS. An example of such a 
calibration plot on the overlay appears in Figure 13. These procedures 
were performed regularly during analyses to insure continual calibration 
of the spectral analyzer and X-Y plotter. 

For purposes of the chart format the RMS voltages were con- 
verted to LOG 10 units and a graduated scale was constructed so that 
the logarithm of volts RMS could be interpolated from spectral plots. 
The amplitude scale was adjusted for each spectrum as a function oF 
both input gain and spectral gain. These values were then converted to 
Popetevels for use in calculating ene values. 

The equation (based on Federal Scientific Specinlcdin taney, 
memcenvert volts RMS to PSD is 

A 


0.79 (cal. level vee? 


(27) 
Filter Bandwidth 


PSD level (v7 /Hz) = 


43 








yy LUN LLL 


Wy 


Mh 


1 y Ni 
VT Wi] 
y, 


ef 
— Wi A we how /| = 


7 Uv 


7an6 






TM, Hy, DAL y) 


ora Namen ete ey 


40 60 cy 


4 


2 


— 
Zz 


6 BR 10 
rz 


iS Cali praelon 


Figure 





where the filter bandwidth (f = 0.4) is a function of the analysis 
range selected. 
b. Obtaining Turbulence Parameters from Scaled Spectra 

The turbulence parameter oe was obtained from the tempera- 
ture variance spectra on the basis of the universal formula, Equation 
(2), for the inertial subrange in wave number space S(k). This expres- 
sion predicts a -5/3 slope for the spectra when plotted in log-log 
formats. 

Figures 14 and 15 are typical spectra considered in the 
analyses. Temperature spectra often exhibited slopes slightly different 
than the expected -5/3. This feature of temperature spectra have been 
observed by others. The existence of "cold spikes" in temperature 
traces has been given as a reason for this, e.g., Friehe (1976). A 
universally accepted reason is still not available, however. 

Assuming -5/3 slopes for the variance spectra, S(f), obtained 
from the analyzer, the intercept of the ~5/3 slope with the l Hz fre- 
quency line was the spectral density denoted (PSD) value used in comput- 
ang the parameter of interest, ales 

The measured PSD value obtained from temperature spectra 


was converted to a spectral density in engineering units by the equation 


Z 
° PSD 
oe 


Sn (f) 


eevee ° Hole RG = 00° Hz (28) 


where or is the calibration factor for the temperature system. Determina- 
tion of the value for C., has been described by Lund (1975). 

Since the temperature fluctuations were measured at a fixed 
point in the flow, the resultant spectral values are defined at 


"temporal" frequencies, £, and are denoted as S, (£) in equation (28). 
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2 
UO Obtain Ch , temporal (f) and space (k) scales had to be 
related in order to use Equation (5). This was accomplished by using 
Taylor's "frozen turbulence" hypothesis discussed in section II,B. The 


following equation relates temporal to wavenumber spectral values, 


7 - 2 ,-2/3 
ig) G8) = ONS Us) ene (29) 


e e e e Z s 
where 5, (f) i the Spectbad aensity value with units of °C /HZ2esihis 
Z ; 
leads to the following relationship between Cr and S(f), which is 


measured, 


& = £S,,(f) (30) 


where, empirically, Cc. =—"0.25 ana K = 2nf/U. 


From Equation (30), with measured values of f, Spf). ane 


OF oe can be determined at each level of interest. These were the 


; bec hak , 2 
procedures used to determine the variation with height of Ch. values, 
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tive RESULTS 


A comparison of the observed results and the Wyngaard et al predic- 
tion curve based on Equation (11) appears in Figure 16 where individual 
data points appear as dots, and averages over Ri intervals of 0.25 
appear as dots within a larger circle. The error bars are standard 
deviations from the mean within each interval, while the number at the 
top of the error bars is the number of obsServations defining the mean 
Value. 

For both the stable (+ Ri) and unstable (- Ri) Stratification cases, 
there appears to be little agreement and no noticeable trends. The 
negative Ri numbers were within acceptable ranges from 0.0 to -2.6. The 
poSitive values, however, ranged from 0.0 to as high as 35.7. Since 
turbulence above a critical Ri number is virtually non-existent due to 
the extreme stability, it was acceptable to assign lower positive Ri 
Values to Ri numbers above 1.0. A critical Ri value is generally con- 
Sidered to be 0.21. 

The scatter in the observed results in Figure 16 can be attributed 
to scatter in both the measured oe values as well as 06/0Z values. 
Deviation of temperature spectra from -5/3 slope caused uncertainty in 
cap estimates. The 06/9Z values are important because the temperature 
gradients are often very small and the quantity is squared for the nor- 
malization. 

In Figure 17, the results computed by Friehe's bulk aerodynamic 
expression, Equation (25), are compared to the Wyngaard et al prediction 


curve. There seems to be some agreement between the observed results 
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Figure 17. Results of the spectrally derived temperature-structure 
parameter vS. A. 
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and the predicted curve close to near-neutral stability conditions, but 
with the large amount of scatter present, this agreement may be coinci-~ 


dental. 


The values obtained for Figure 16 and Figure 17 are summarized in 


Mewes Li and Iil. 
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Ve. CONCLUSIONS 


Correlation of spectrally derived ce results with the stability 
parameter zZ/lL or with the Richardson number values was not good for 
these data in a marine environment. The temperature gradient ina 
marine environment is too small and gives rise to very large om values 
which contribute to the observed scatter. As pointed out by Friehe 
(1976), a complication in the marine boundary layer is the sometimes 
anomalous behavior of the temperature field where inertial subranges 
are not always obtained. Salt spray is also a factor to content with 
in a marine environment and can be the cause of many problems with 
instruments exposed to it. 

For Navy operational purposes, then, existing empirically derived 
expressions defined by z/L or the Richardson number may require modi- 
fication before use in describing the atmosphere over an ocean environ-~ 
ment for obtaining coe indirectly. However, for reasons explained in 
the preceding paragraph, more measurements using a different approach 
iecomputing se may lead to better correlation between ca and an 
easily defined stability parameter. It should also be noted that the 
correct empirical relationship, when finally derived, may be different 
over the marine environment than the curve predicted by Wyngaard et al 


over land. 
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